ApJ Letters, in press 

Preprint typeset using IATf^X style cmulateapj v. 14/09/00 



THE EVOLUTION OF FIELD EARLY-TYPE GALAXIES TO Z ~ 0.7 

1 

TOMMASO TREU 

California Institute of Technology, Astronomy Department, 105-24, Pasadena, CA 91125 

tt@astro.caltech.edu 
AND 

Massimo Stiavelli, S. Casertano 

Space Telescope Science Institute, 3700 San Martin Dr., Baltimore, MD 21218 

AND 

PALLE M0LLER 

European Southern Observatories, Karl-Schwarzschild Str. 2, D85748, Garching bei Miinchen, Germany 
1> ' AND 

o : 

►7- , Giuseppe Bertin 

Universita degli Studi di Milano, Dipartimento di Fisica, Via Celoria 16, 120133, Milano, Italy 

' ApJ Letters, in press 



(N 



O 



6 



ABSTRACT 



We have measured the Fundamental Plane (FP) parameters for a sample of 30 field early- type galaxies 
(E/S0) in the redshift range 0.1 < z < 0.66. We find that: i) the FP is defined and tight out to the 

bio-boat rorlcbift bin- ii ) tbo !n+prppnt *v pvnluoc nc ri^j I rl z> — .^^^"9"^? 



highest redshift bin; ii) the intercept 7 evolves as d-y/dz — 0.58_q '-^3 (for f2 = 0.3, J7a = 0.7), or, in terms 



, of average effective mass to light ratio, as d\og(M/LB)/dz = — 0.721°'^, i. e. faster than is observed for 

cluster E/S0 (—0.49 ± 0.05). In addition, we detect [Oil] emission > 5A in 22% of an enlarged sample 
of 42 massive E/S0 in the range 0.1 < z < 0.73, in contrast with the quiescent population observed in 
clusters at similar z. We interpret these findings as evidence that a significant fraction of massive field 
E/S0 experiences secondary episodes of star-formation at z < 1. 



Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution — - galaxies: formation 
— galaxies: structure 

1. introduction both in clusters and in the field. In addition, hierarchi- 



The Fundamental Plane (Djorgovski & Davis 1987; cal clustering models (e g Kauffmann 1M6) predict the 

. . , -p. , , , ,„„, , r, • i- 1.4. 1 a- u stellar populations of held E/ SO to be significantly younger 

> Dressier et al. 1987; hereafter IP) is a tight correlation be- ,. „ , , „ „, „ ° ,.~ ° 

t, i.u a i.- j- /D wv a i° c u • 1.4. than the ones of cluster E/S0. 1 he effects of age differences 

• i-H , tween the effective radius (R e ), the effective surface bright- , . . ' , , TT . & ,„ .. , 

k>( /en \ i ii , \ ' , . . / \ are difficult to observe m the local Universe Bcrnardi et 

X] ■ ness (Sli e ), and the central velocity dispersion (a) , ..„„> . , . . . v , . 

al. 1998), when the typical population ages are large, but 

03 logi? e = alog fi + SB e + 7 (1) &re greatly enhanced at intermediate redshift. For these 

that is observed to hold for local early-type galaxies reasons, we have embarked on a campaign aimed at mea- 

(E/S0). Under the assumption that E/S0 are homologous suring the evolution of the FP of field E/S0 with redshift. 

dynamical systems the FP can be interpreted in terms of i n previous papers (Treu et al. 1999; 2001a, b; hereafter 

a power law relation between mass (M) and mass-to-light T01a,b), we have analyzed a sample of 19 E/S0, finding 

ratio (M/L) (Treu et al. 2001b and references therein). In that a tight FP exists in the field out to z « 0.4 and that 

recent years, it has been found that a tight FP exists in the evolution of the intercept is marginally faster than in 

clusters out to redshift z = 0.83 (van Dokkum & Franx the clusters, indicating a marginally younger age for field 

1996; Pahre 1998; Bender et al. 1998; Kelson et al. 2000). galaxies (see also van Dokkum et al. 2001 and Kochanek 

The modest evolution of the intercept 7 and the tightness e t al. 2000). 

of the FP at O.f < z < 1 can be explained in terms of Here we present results from a larger sample of 30 E/S0 

passive evolution and an old age of the stellar populations extended to z ~ 0.7, increasing dramatically our sensitiv- 

in cluster E/S0. The absence of a dramatic evolution in ity to differences between cluster and field environment, 

the slopes a and (3 argues against the interpretation of the anc i s h w that field E/S0 evolve faster than cluster ones 

FP "tilt" resulting solely from a mass-age relation. at the 95% CL. In addition, we report on the detection 

So far, most of the studies have been focused on the clus- f [OII]3727 emission in a significant fraction of massive 

ter environment. However, the population of galaxies in ne ld E/S0 at intermediate redshift, which we interpret as 

clusters is likely to evolve with redshift by accretion of field evidence for secondary episodes of star- formation, 

galaxies. Therefore, in order to obtain a complete and reli- We assume the Hubble constant, the matter density, 
able picture, it is necessary to study the evolution of E/S0 

1 Based on observations collected at ESO (Paranal) under programmes 65.O-0446 and 66.A-0362, and with the NASA/ESA HST, obtained at 
STScI, which is operated by AURA, under NASA contract NAS5-26555. 
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and the cosmological constant to be respectively Ho = 
50/i5okm s _1 Mpc _1 (/i 50 = 1.3 when necessary), = 0.3, 
and fl A = 0.7. 

2. SAMPLE SELECTION AND OBSERVATIONS 

The sample of E/SO was selected on the basis of 
Hubble Space Telescope (HST) images taken from the 
HST-Medium Deep Survey (MDS; Griffiths et al. 1994). 
The selection criteria are identical to the ones ex- 
tensively discussed in TOla, except for the color cut 
(1.25 < V606-I814) chosen to select higher redshift E/SO; 
V606 and I§i4 indicate Vega magnitudes through HST fil- 
ters F606W and F814W respectively) and the magnitude 
range (19.3 < Igi4 < 20.3). The effects of color and mag- 
nitude selection are taken into account in the analysis pre- 
sented here as discussed in TOlb. Images are available at 
the MDS web site at http://www.archive.stsci.edu/mds. 
Morphological classification was performed independently 
by two of us (TT, MS), based on visual inspection of the 
images, of the residuals from the fit, and of the luminosity 
profiles. In particular, since contamination from Sa galax- 
ies 2 could bias our results towards a younger age of the 
stellar populations, we were extremely conservative in re- 
jecting galaxies that showed sign of a spiral disk either in 
the direct images (F606W and F814W), or in the residuals 
from the r 1 / 4 fit. As independent check of our classifica- 
tion, the bulge-to-total luminosity ratio measured by the 
MDS group for our E/SO is always larger than 0.6 and 
typically in the range 0.8-1. 



Ol=ust00 15 7=0.646 <j=200±18 Pl=ust00 16 z=0.554 (t=145±18 




A(A— rest frame) 



Fig. 



1. — Spectra of galaxies with measured central velocity dis- 
persion. The shaded regions are affected by sky-line residuals and 
atmospheric absorption and have been masked out during the fit. 
Note [Oil] in emission in galaxies B2 and A5. 



a resolution of ~ 90 — 100km s ; exposure times ranged 
between 2 x 1800s and 2 x 3600s. For 10 out of 16 E/SO 
observed (Figure 1) we obtained a reliable central velocity 
dispersion (the success rate depending on observing con- 
ditions and redshift, since the observed wavelength range 
was optimized for the range z « 0.45 — 0.65). The data 
reduction was very similar to the one described in TOla. 
Combined with the data of TOlb, our total sample consists 
of 42 galaxies (in the range z = 0.10 — 0.73), 30 of which 
with measured velocity dispersion (z = 0.10 — 0.66). 

3. THE EVOLUTION OF THE FUNDAMENTAL PLANE 

Panels (a) to (e) in Figure 2 show the location in the 
FP-space of the galaxies in our sample, binned in redshift. 
The data (open symbols) in the rest-frame B band 3 are 
compared to the local relation shown as solid line. The 
main result is that the brightening at fixed a and i? c in- 
creases with redshift. By assuming constant slopes a and 
(3 we obtain the average evolution of the intercept with 
redshift plotted as filled pentagons in panel (f). 
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Spectra for 16 galaxies were obtained in service mode 
from April 2000 to March 2001, using the Focal Reducer 
and Spectrograph 2 (FORS2) at the Very Large Telescope 
(VLT). The adopted grism R600 with a 1" wide slit gives 

2 See, e.g., Smail et al. (1997) and van Dokkum et al. (1998a) for a discussion. 

3 The analysis in the V band, or using different Coma FP, yields the same results, see TOlb for discussion 



Fig. 2. — FP in the rest-frame B band. In panels (a) to (c) we 
show the field E/SO (open stars if [Oil] emission is detected, open 
pentagons otherwise), binned in redshift and compared to the FP 
found in the Coma Cluster (Bender et al. 1998). Panel (f) shows 
the average offset of the intercept of field galaxies from the local 
FP relation as a function of redshift (large filled pentagons), com- 
pared to the offset observed in clusters (open squares; van Dokkum 
& Frame 1996; Kelson et al. 1997; Bender et al. 1998; van Dokkum 
et al. 1998; Kelson et al. 2000). The solid lines represent the evo- 
lution predicted for passively evolving stellar populations formed in 
a single burst at z = 1,2,5 (from top to bottom) computed using 
Bruzual & Chariot (1993) models in the BC96 version (as described 
in TOlb). We assume U = 0.3, U A = 0.7 and h s0 = 1.3 

Under the assumption of passive evolution and constant 
slopes (TOlb), the evolution of 7 can be converted into 
the evolution of the average M/L ratio, which can be read 
on the right scale of panel (f). Modelling the evolution as 
"f(z) = 7(0) +7'(0)z, where the prime indicates derivative 
with respect to z, a least \ 2 fit yields 7^ = 0.64. How- 
ever, in order to avoid biased results, it is crucial to take 
into account the selection process, for example by apply- 
ing the Montecarlo-Bayesian method introduced in TOlb, 
generalized by allowing the scatter to vary as a function 



The Fundamental Plane of field early-type galaxies 



of redshift, cr-y(z) — 07 (0) + cr'(0)z. In this way, we derive 
the posterior probability density given the set of obser- 
vations, p(7', a' |{7i}) in the notation of TOlb, shown in 
Figure 3 panel (a). The posterior probability peaks at 
i B = 0.58 (0.45 - 0.67 68% limit), which corresponds to 
\og{M/L B )' = —0.72+2^6- Tne scatter is constant or at 
most mildly increasing with redshift. Note that the effects 
of the magnitude selection limit (see TOlb) on the mea- 
sured evolution of the intercept depend on the intrinsic 
scatter of the FP, and therefore it is important to study the 
combined probability density. Note also that the galax- 
ies with Oil emission (open stars in panels c and d) fall 
slightly above the quiescent ones (pentagons), as expected 
if their M/L is smaller because of a fraction of young stars. 
However, they do not dominate the general trend, since the 
results change very little if they are removed from the sam- 
ple (e.g. 7^ from the least \ 2 fit goes from 0.64 to 0.61). 
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Fig. 3. — Panel (a): probability contours for the variation of the 
intercept and scatter with redshift (respectively 7' and ai^ ; the pen- 
tagon marks the peak of the probability density, the dashed line the 
68% limit and the solid line the 95% limit). Panel (b): probabil- 
ity contours for the variation of the intercept (7') and slope (a') 
with redshift, obtained by assuming a' = 10/3'; symbols as in panel 
(a). Panel (c): evolution of M/Lg as a function of M (denned as 
5cr 2 R e /G), as obtained by assuming fixed slopes (dashed line), or 
the values of a' and f3' corresponding to the peak of the probability 
density in panel (b; solid line). 

It is desirable to generalize the treatment by allowing 
also the slopes to evolve with redshift. Even though a 
larger sample is needed to this aim, possibly spanning a 
wider range of parameters, we will try to address this is- 
sue by assuming two additional constraints: i) the scatter 
does not change from z = to z ~ 0.7 and ii) a' = 10/3' 
(hence the FP is equivalent to a relation between M and 
M/L at any redshift). Under these assumptions, we ob- 
tain the results shown in Fig. 3 panels (b) and (c), plotted 
using M = 5cr 2 R e /G (Bender, Burstein & Faber 1992). 
The evolution of the slopes remains highly unconstrained 
(panel b), but the evolution of the M vs M/L relation is 
better constrained (panel c). For the typical masses in our 
sample (< log M/M >= 11.4) the derived value of M/L' 
does not depend significantly on the value of a', provid- 
ing a test of the robustness of the measurement. On the 
other hand, the result shown as a solid line in panel (c) 
would indicate that the evolution of M/L with redshift 
becomes slower as the mass increases (i.e. for example 
massive E/S0 are older). More data are needed to confirm 
these findings. 

In Figure 2 we also show the evolution of 7 predicted 
by assuming passively evolving single stellar populations. 
The best-fitting single burst redshift of formation should 
be interpreted in terms of a luminosity weighted redshift 
of formation, since it is clear that a more complex scenario 
than passive evolution of single coeval stellar populations 
is needed to explain the present observations (TOlb). Fi- 
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nally, cluster data taken from the literature are shown as 
open squares in Figure 2 for comparison. The field FP 
evolves faster than the cluster one: for the cluster van 
Dokkum et al. (2001) measure \og{M/L B )' = -0.49±0.05, 
i. e. 7' = 0.39 ± 0.04, which falls on the 95% contour in 
panel (a) of Figure 3. 

4. EMISSION LINE PROPERTIES 

At z > 0.3 the emission line [OII]3727 is redshiftcd into 
the wavelength range covered by our instrumental setup. 
Six out of 27 E/S0 for which [Oil] is detectable show sig- 
nificant emission (rest-frame equivalent width EW> 5 A). 
If attributed to star-formation, the observed [Oil] EW 
correspond to star- formation rates of order 0.5-5 M^/yr 
(Kennicutt 1992). The fraction of E/S0 with sizable star- 
formation in our sample is significantly larger than the one 
found for massive E/S0 in clusters at similar redshifts and 
in the local Universe. In Figure 4 we plot the fraction 
of E/S0 with [Oil] EW> 5 A as a function of redshift for 
field (pentagons) and cluster (squares) samples taken from 
the literature. 
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FlG. 4. — Fraction of E/S0 with [Oil] emission larger than 5 
A EW (rest-frame) as a function of redshift. Pentagons repre- 
sent field galaxies, squares cluster galaxies. The points arc filled 
if selected similarly to our sample. Field galaxies are from Cald- 
well (1984), Kennicutt (1992), Janscn et al. (2000; limited to 
Wig— 5 log /150 < —20, corresponding roughly to our survey), Brinch- 
mann et al. (1998; limited to lg < 20.3), and the present work (large 
filled pentagon with error bars). The filled squares arc from Fisher 
et al. (1998), van Dokkum et al. (1998a), van Dokkum et al. (2000). 
The open squares are from Smail et al. (1997) and Dressier et al. 
(1999). 

E/S0 with similar [Oil] EW, or similarly with relatively 
blue colors, have been reported previously, suggesting that 
secondary episodes of star-formation are common at in- 
termediate redshift (Schadc et al. 1999; Mcnantcau et al. 
2001). Im et al. (2001) measured the width of the [Oil] 
line for a sample of blue E/S0 at intermediate redshift 
finding velocity dispersions a < 80 kms" 1 ; based on this 
measurement they argue that most blue E/S0 are not the 
progenitors of present-day massive E/S0, but rather less 
massive spheroids, for which star-formation is observed 
also in the local Universe (Jansen et al. 2000). Our sam- 
ple adds further information, because it is made of bright 
(I < 20.3) and massive E/S0. The velocity dispersions 
we obtained via absorption-line kinematics are typical of 
massive field E/S0. In particular, the three objects with 
significant [Oil] emission for which velocity dispersion is 
available have a = 118, 179, 233 kms" 1 . From the FP rela- 
tionship, as measured from our sample, we estimate for the 
other three galaxies with [Oil] emission a = 105, 171,261 
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kms -1 . 

We can estimate the amount of stellar mass assembled 
in these secondary burst in the following way. The proba- 
bility of observing a burst is given by 

p(OII)=n^, (2) 

where n is the average number of bursts per galaxy i is 
the average duration of the burst, and At is the interval 
in cosmic time between z — 0.73 and z — 0.26 when the 
burst was observable. The total average mass in secondary 
bursts is then: 

M on = ntM = p(Oll)AtM. (3) 
By using the observed values for p(OII) = 6/27 and 

M — 0.5 — 5M Q yr _1 , we estimate that the average stellar 
mass formed in these bursts is of the order of Mqjt ~ 
5 • 10 8 - 5 • 1O 9 M . 

5. CONCLUSIONS 

Our sample of 30 field E/SO defines a tight FP out to 
z = 0.66, with no or modest increase of the scatter with 
rcdshift. The intercept 7 evolves with redshift as 7' = 0.58, 
which, interpreted in terms of passive evolution of the 
stellar populations, implies log(M/Ls)' = — 0.72+g 16 , i.e. 
faster than that observed in clusters by other groups (e.g. 
-0.49±0.05 van Dokkum ct al. 2001). In addition, 22% of 



the galaxies show [Oil] in emission, with EW > 5A. As- 
suming the emission is due to star-formation, we estimate 
that such bursts produce of order Mqjj <~ 5-10 8 — 5-10 9 M Q 
of stellar mass between z = 0.73 and z = 0.26. 

The existence and tightness of the FP suggest that no 
major structural changes occur between z ~ 0.7 and z = 0. 
However, we find evidence that some stellar mass is formed 
at relatively recent time during secondary bursts. Al- 
though these bursts contribute only a small fraction of 
the total stellar mass they contribute significantly to the 
evolution of the observable luminous component. For ex- 
ample, a scenario where most of the stars formed at z > 1 
and secondary star-formation occurs at z < 1 could not 
only explain the evolution of the FP (TOlb) and the ob- 
served [Oil] emission but could also reconcile the mod- 
est evolution of the number density of E/S0 observed to 
z ~ 1 (Schade et al. 1999; Im et al. 2001, in press), with 
the paucity of red E/S0 observed in some infrared surveys 
(Menanteau et al. 1999; Treu & Stiavelli 1999; see also 
Jimenez et al. 1999 and Daddi et al. 2000). 
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Fitting Software developed by R. van der Marel and 
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